Mon. Not. R. Astron. Soc. 000, 000-000 (2009) Printed 19 April 2010 (MN MfeX style file v2.2) 



Cluster Galaxies Die Hard 



Simone M. Weinmann 1 *, Guinevere Kauffmann 1 , Anja von der Linden 2 , 
Gabriella De Lucia 3 

1 Max-Planck Institut fuer Astrophysik, Karl-Schwarzschild-Str.l, Postfach 1317, 85741 Garching, Germany 

2 Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA 94305-4085, USA 
3 Osservatorio Astronomico di Trieste INAF, Via Tiepolo 11, 34143 Trieste, Italy 



o 

(N 

Or 
< 



< 



ABSTRACT 

We investigate how the specific star formation rates of galaxies of different masses 
depend on cluster-centric radius and on the central/satellite dichotomy in both field 
and cluster environments. Recent data from a variety of sources, including the cluster 
catalogue of von der Linden et al. are compared to the semi-analytic models of De 
Lucia & Blaizot. We find that these models predict too many passive satellite galaxies 
in clusters, too few passive central galaxies with low stellar masses, and too many 
passive central galaxies with high masses. We then outline a series of modifications 

■ to the model necessary to solve these problems: a) Instead of instantaneous stripping 
Q-T of the external gas reservoir after a galaxy becomes a satellite, the gas supply is 
q assumed to decrease at the same rate that the surrounding halo loses mass due to tidal 
5-h ■ stripping, b) The AGN feedback efficiency is lowered to bring the fraction of massive 
£jj ' passive centrals in better agreement with the data. We also allow for radio mode AGN 

feedback in satellite galaxies, c) We assume that satellite galaxies residing in host 
haloes with masses below 10 12 /i _1 Af© do not undergo any stripping. We highlight 
the fact that in low mass galaxies, the external reservoir is composed primarily of 
gas that has been expelled from the galactic disk by supernovae driven winds. This 
gas must remain available as a future reservoir for star formation, even in satellite 
galaxies. Finally, we present a simple recipe for the stripping of gas and dark matter 

■ in satellites that can be used in models where subhalo evolution is not followed in 
CN \ detail. 
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1 INTRODUCTION feedback, feedback by active galactic nuclei, and morpholog- 
ic ■ ical transformations. All of these processes are likely to be 
^ It has long been known that cluster galaxies are redder, less . a , , , r ,. , . 
jrt ' influenced by the decline ot gas accretion in galaxies residing 
active in their star formation, and of earlier type than galax- . , , . ™ , , „ 

„ , , , _ , ' „ , ;„„„ „ , ° m group and cluster environments. Clusters and groups ot 

ies in the field (e.g. Oemler 1974; Dressier 1980; Balogh et 



al. 1997). This difference has been found to become more 



galaxies thus provide unique laboratories in which we can 
observe the evolution of galaxies under conditions different 
pronounced towards the centers of clusters (e.g. Postman ^ those in the field 
& Geller 1984; Goto et al. 2004; Barkhouse et al. 2009; ° ° Se m e e ' 

Hansen et al. 2009; von der Linden et al. 2009) and per- Recently, several studies have investigated the fraction 

, , . ! . / m . . of passive cluster galaxies as a function of cluster-centric 

haps also with increasing cluster mass (e.g. Wemmann et v ° 

al. 2006a; Martinez et al. 2006, Kimm et al. 2009; but see radius at fi xed steUar mass - Most of them find that an in " 

also De Propris et al. 2004 and Tanaka et al. 2004). Under- crease of the P assi , ve fraction towards the center of cluster 
, u ■ , , rr , „ • , , • , c is still detectable (Bamford et al. 2008; von der Linden et 

standing these environmental effects is not only crucial lor v ' 

modelling the detailed properties of the global galaxy pop- aL 2009 ^ this work ^ but see also van den Bosch et aL 2008 )' 
, . • i , ., i u i i • i . s. c j .1 Investigating the trends at fixed stellar mass is important 

ulation, but it can also hefp us to investigate fundamentaf e> & r 

c , , ,. ,., , c -t- since the galaxy stellar mass function may depend on en- 

processes of galaxy evolution like star formation, supernova , 

vironment (e.g. Balogh et al. 2001; Baldry et al. 2006; but 

see also von der Linden et al. 2009), and since galaxy prop- 
erties are strongly correlated with stellar mass (Kauffmann 
* E-mail:simone@MPA-Garching.MPG.DE et al. 2003). In this work, we determine the same relations 
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in the cluster catalogue of von der Linden et al. (2007), in 
which the brightest cluster galaxies are selected with great 
care and are thus likely to mark the approximate center of 
the cluster. This is crucial for deriving correct trends as a 
function of cluster-centric radius. 

Semi-analytical models (hereafter referred to as SAMs, 
e.g. Kauffmann et al. 1993; Cole et al. 2000; Croton et al. 
2006) track the evolution of galaxies over time, using sim- 
ple prescriptions for the treatment of gas-physical processes 
combined with analytical merger trees of dark matter haloes 
or trees derived from N-body simulations. Most of these 
models treat environmental effects in a simplistic manner. 
Motivated by the "starvation" scenario suggested by Larson, 
Tinsley & Caldwell (1980), all hot gas around the satellites 
is immediately removed upon infall. The stripped gas is then 
made available for cooling to the central galaxy of the corre- 
sponding Friends-of-Friends group. This simple prescription 
leads to satellite galaxies that are too red (e.g. Weinmann 
et al. 2006b, Wang et al. 2007) compared to observations. 
This suggests that part of the hot gas should remain with 
the satellite galaxy. This is not a minor issue, since satellite 
galaxies constitute a significant fraction of the total galaxy 
population and since changing the prescription for gas strip- 
ping in satellites may have a considerable impact on the 
central galaxy population. First, satellites eventually merge 
with the central galaxies in their halo. If they can grow to 
higher stellar masses, central galaxies will become more mas- 
sive as well. Second, if part of the hot gas stays attached to 
satellites, the amount of gas available for cooling to the cen- 
tral galaxy is reduced. Third, if satellite galaxies merging 
with the central galaxy still contain cold gas, the ensuing 
star burst will make the central galaxy bluer for a certain 
period of time, and will result in a higher final stellar mass. 

Attempts to treat environmental effects in SAMs more 
realistically have already been made by Kang & van den 
Bosch (2008) on the basis of the Kang, Jing & Silk (2006) 
SAM and by Font et al. (2008), on the basis of the Bower 
et al. (2006) model. Kang & van den Bosch (2008) found 
that their simple prescription for decreasing the efficiency 
of gas stripping in satellites leads to a fraction of blue cen- 
tral galaxies which is higher than observed. They suggested 
counterbalancing this effect by the inclusion of an additional 
prescription for the disruption of satellites. Font et al. (2008) 
implemented a more sophisticated model based on the hy- 
drodynamical simulations of ram-pressure stripping by Mc- 
Carthy et al. (2008), and obtained better agreement with 
observed environmental effects than previous semi-analytic 
models. 

Ram-pressure stripping of the cold disk gas in satellites 
has been studied in the SAMs of Okamoto & Nagashima 
(2003) and Lanzoni et al. (2005). These studies concluded 
that this effect has a negligible impact on the results, be- 
cause the complete stripping of the hot gas halo already 
makes the satellites passive. 

In this work, we modify the SAM of De Lucia & Blaizot 
(2007, hereafter DLB07) in order to reproduce (i) the rela- 
tion between passive fraction and cluster-centric radius, as 
well as the passive fraction in field galaxies and (ii) the dis- 
tribution of specific star formation rates of satellite galax- 
ies. Up to now, this combination of observational relations 
has never been used to constrain SAMs (but see Diaferio 
et al. 2001, who study passive fractions as a function of 



cluster-centric radius in a set of SAMs implemented on low- 
resolution N-body simulations) . We apply a realistic cluster 
finder to the DLB07 SAM to allow for a fair comparison with 
observations. We find that a very simple model in which the 
diffuse gas halo around satellites is stripped at the same rate 
as the dark matter subhalo loses mass due to tidal effects 
gives good agreement with observations. Interestingly, we 
find that this stripping does not proceed exponentially, but 
linearly, and that this behaviour is crucial for reproducing 
the distribution of the specific star formation rates (SSFR) 
in satellites. We also implement a series of modifications to 
the SAM, which lead to improved agreement with the ob- 
served passive fractions as a function of stellar mass for the 
central galaxies. Our goal is not to obtain a full model which 
reproduces all the observations; rather, we look at a variety 
of models which allow us to understand better the com- 
plex interplay between feedback by supernovae ('SN' here- 
after), feedback by active galactic nuclei ('AGN' hereafter) 
and satellite galaxy stripping. We also test a model in which 
the hot gas of satellite galaxies is removed by ram-pressure 
stripping, but find that this model does not reproduce ob- 
servations in detail. 

To summarize, the main goal of this study is to pro- 
vide insight into what causes environmental effects, on which 
timescales they act and how they might be modelled in 
SAMs of galaxy formation. In section [2] we present the ob- 
servational and semi-analytical data and explain the con- 
struction of our mock cluster catalogue. In section [3] we 
present the observational results which we use to compare 
to models. In section[4] we discuss the stripping of dark mat- 
ter subhaloes. In section [5] we present modified versions of 
the SAM of DLB07, which we compare with observations in 
section [5] Finally, in section [7J we discuss our findings and 
test a model of ram-pressure stripping for comparison. In 
section [8j we give a summary of our results. 



2 DATA 

2.1 The Cluster Catalogue 

We use the Cluster Catalogue of von der Linden et al. (2007, 
hereafter vdL07) which is based on the SDSS DR4 and the 
C4 Cluster Catalogue (Miller et al. 2005). The C4 Cluster 
Catalogue identifies clusters in a parameter space of posi- 
tion, redshift and colour, making use of the fact that at 
least a fraction of the cluster galaxies lie on a tight colour- 
magnitude relation. VdL07 carefully identified the brightest 
cluster galaxies (BCGs) in these clusters, making sure that 
galaxies that were not targeted spectroscopically would not 
be missed. They then redetermined cluster memberships and 
velocity dispersions using the bi-weight estimator by Beers 
et al. (1990), and estimated cluster masses from the velocity 
dispersions, as described in vdL07. Their final sample con- 
sists of 625 clusters at redshifts between 0.03 and 0.1, with 
masses between 10 12 and 10 15 Mq. In most of what fol- 
lows, we will focus on clusters with masses 
of which there are 214 in the sample. The 341 clusters with 
masses between 10 13 — 10 14 h^ 1 Mq will only be used for com- 
parison. At the redshift of the most distant cluster, galaxies 
were observed down to a limiting magnitude of M r - 5\ogh— 
-19.75. 
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In all of what follows, we weight each galaxy by the 
inverse of the maximum volume out to which it can be ob- 
served, given the apparent magnitude limit of the survey. 
Stellar masses are calculated using the method of Kauff- 
mann et al. (2003). 

2.2 Central Galaxies 

In order to obtain large observational samples of central 
galaxies needed for comparison to the SAM, we make use 
of the SDSS DR4 group catalogue of Yang et al. (2007). 
This group catalogue has been constructed by applying 
the halo-based group finder of Yang et al. (2005) to the 
New York Value-Added Galaxy Catalogue (NYU-VAGC; 
see Blanton et al. 2005). From this catalogue, Yang et al. 
(2007) selected all galaxies in the Main Galaxy Sample with 
an extinction corrected apparent magnitude brighter than 
m r — 18, with redshifts in the range 0.01 < z < 0.20 
and with a redshift completeness C z > 0.7. The group 
catalogue is publically available and can be downloaded 
from http: //www. astro .umass . edu/~xhyang/Group.html. 
We refer the reader to Yang et al. (2007) for a more detailed 
description. We only use central galaxies from this catalogue 
that have redshifts in the range 0.01 < z < 0.1 in this work. 

2.3 The simulation and the SAM 

We use the SAM by DLB07 which is based on the Mil- 
lennium Simulation (Springel et al. 2005). The SAM uses 
analytical prescriptions for gas accretion, gas cooling, star 
formation, SN feedback, AGN feedback, dynamical friction, 
merger- induced star bursts and reionization. These are im- 
plemented on the merger trees extracted from the Millen- 
nium Simulation, which follow dark matter haloes and sub- 
haloes over time. Photometric properties of galaxies are 
computed using models for stellar population synthesis and 
dust. More details can be found in DLB07, Croton et al. 
(2006) and references therein. The SAM has been tuned to 
reproduce key observations like the luminosity function at 
z=0 (Croton et al. 2006; DLB07). The aspect of the SAM on 
which we will mainly focus in this paper is the treatment of 
gas in satellite galaxies, which will be explained in more de- 
tail in section [STT1 The cosmological parameters used in this 
model are in agreement with the WMAP1 data (Spergel et 
al. 2003), with U M = 0.25, Q A = 0.75, fi 6 =0.045 and <t 8 =0.9. 

2.4 The SAM Cluster Catalogue 

We now describe the construction of our mock cluster cata- 
logue based on the 2=0 output of the Millennium Simulation 
combined with the DLB07 SAM. Our goal is to mimic the 
procedure used by vdL07 to determine cluster memberships 
and velocity dispersions. However, we do not try to mimic 
the initial identification of clusters in the C4 Cluster Cat- 
alogue, or the method for idenitifying the brightest cluster 
galaxies (BCG). We make two assumptions: first, that if 
we impose the same magnitude limit on the cluster galaxy 
sample from the Millennium Run simulation as we do for 
the vdL07 catalogue, the resulting sample of simulated clus- 
ters can be compared directly with clusters of the same mass 
from our observational sample. Second, we assume that the 



BCGs in vdL07 are correctly identified. The first assump- 
tion is likely not quite correct; our sample probably contains 
more low mass clusters than the sample of vdL07, since those 
are more difficult to identify observationally. The second as- 
sumption seems reasonable. VdL07, Best et al. (2007) and 
Koester et al. (2007) show that the BCGs of vdL07 differ 
systematically from other galaxies of the same mass and are 
good tracers of the cluster centers defined using X-ray im- 
ages of the hot intracluster gas. 

Clusters in the Millennium Simulation are identified 
with a Friends-of-Friends (fof ) algorithm (Springel et al. 
2005). These fof clusters are not necessarily spherical, but 
are often elongated in shape. We select the FOF clusters with 
10 13 < M 2O o,Mm/(/i" 1 M ) < 10 15 , with M 200 ,Miii the mass 
enclosed within a sphere with a density of 200 times the crit- 
ical density of the universe. Subhaloes are identified using 
the SUBFIND algorithm (Springel et al. 2001). By definition, 
the "central galaxies" of clusters reside at the center of mass 
of the most massive subhalo in the FOF group. We use these 
"central galaxies" as starting points for determining cluster 
memberships and velocity dispersion, as done in vdL07 for 
their observational sample of BCGs. 

We mimic the observational selection of vdL07 by only 
using galaxies with M r — 51og/i < —18.6 for the construction 
of our mock cluster catalogues, which roughly corresponds 
to the absolute magnitude limit at the median redshift of the 
vdL07 sample. We add the Hubble flow to the z-direction 
of the velocities of galaxies in our simulation box, and we 
place all galaxies at a minimum distance of 60 Mpc from the 
virtual observer, so that all clusters are at least as far away 
as the nearest cluster in the sample of vdL07. We then follow 
the procedure explained in detail by vdL07, in which the bi- 
weight estimator of Beers et al. (1990) is applied iteratively 
to galaxies surrounding the cluster center, with the velocity 
dispersion and the cluster redshift being redetermined at 
each iteration step. If a galaxy belongs to more than one 
cluster in the final list, then we follow von der Linden et al. 
(2009) and only count it as a member of the one cluster for 
which the quantity 

£>efi = x/{R/R2oo) 2 + (V7a 1D )2 (1) 
is smaller, with 

R 2 = (x- x cl f + (y - y cl ) 2 (2) 
and 

V 2 = (Weff.cl ~ Uoff.gal) 2 - (3) 

We only deviate from vdL07 in two ways : 

• We define two cluster mass bins, one with M2oo,mhi = 
10 13 - lO 14 ft _:L M and one with M 2 oo, M m = 10 14 - 
10 15 Mq. If a galaxy lies close to the centers of two clus- 
ters in the same mass bin, we only allow it to enter into 
the iterative process for the closer cluster, with the distance 
defined by equation Q] This prevents cluster galaxies from 
being falsely linked to a neighbouring cluster of similar mass. 
Such a step is not included in vdL07, where a given galaxy 
can in principle enter into the iterative process for several 
nearby cluster centers. We believe that this step is justified 
here, as the density of clusters in the simulation is higher 
than in an observational sample. We have checked that this 
constraint only increases the fraction of interlopers in the 
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outskirts of cluster by around 10 % and does not affect our 
final results. 

• The initial guess for the ID velocity dispersion is taken 
to be 500 km/s for all clusters. VdL07 calculate an ini- 
tial guess individually for each cluster, using the average of 
the different velocity measurements in the C4 Cluster Cat- 
alogue. We have checked that using an initial guess of 250 
km/s or 1000 km/s makes virtually no difference to our final 
results. 

To enable faster data processing, we only use a subvol- 
ume of the Millennium Simulation, which encloses around 
1/5 of the total simulation box and contains the centers 
of 7586 clusters with masses 10 13 < M 200 .Mill 
10 15 . We have checked that the mass distribution of these 
clusters is the same as in the entire simulation box. 412 
clusters have a mass above 1O 14 /i _1 M0. After applying our 
bi-weight estimator to all of the 7586 central cluster galax- 
ies as described above, we remove all clusters which contain 
less than 4 members in the final iteration step, as done by 
vdL07. We also remove (i) all clusters which are clearly offset 
from the relation between velocity dispersion and number of 
member galaxies defined by the bulk of the cluster popula- 
tion and (ii) all clusters for which the position of the original 
BCG and the cluster center determined by the bi-weight es- 
timator of Beers et al. (1990) differ by more than 0.002 in 
redshift. This step is necessary since such clusters have un- 
dergone special treatment in vdL07, which we do not repeat 
here. It affects less than 3 % of the clusters in our sample. 
Finally, we are left with 4393 clusters, constituting what will 
be called the "Beers- Millennium Cluster Catalogue" in what 
follows. The reason for the high reduction in sample size 
going from the Millennium Cluster Catalogue to the Beers- 
Millennium Catalogue is simply that our initial Millennium 
Cluster sample contains many relatively low mass clusters, 
of which many contain less than four galaxies. The ID ve- 
locity dispersion of clusters with a low number of member 
galaxies also has a large scatter, which makes convergence 
of the bi-weight estimator of Beers et al. (1990) less likely. 

We estimate M200 for the clusters in our final sam- 
ple from the ID velocity dispersions according to eq. 1 of 
vdL07. In all of what follows, M200 or the "cluster mass" 
refers to the mass as obtained in this way. 3151 of the clus- 
ters in the "Beers-Millennium Cluster Catalogue" are as- 
signed a mass above 10 13 /i _1 M Q , 747 above 10 14 fr _1 M©. 
We therefore have significantly more galaxies with masses 
above 10 14 fo _1 Mq than in our initial Millennium Cluster 
Catalogue. This already indicates that a non- negligible num- 
ber of clusters have masses that are overestimated. 

We now compare the properties of the Millennium Clus- 
ter Catalogue catalogue and the Beers-Millennium Clus- 
ter Catalogue in detail. In the top panel of Fig. Q] we 
show the distribution of halo masses in the three cata- 
logues used in this work: The vdL07 cluster catalogue, 
the original Millennium cluster catalogue and the recov- 
ered Beers-Millennium catalogue. The original distribution 
of halo masses is clearly smeared out for the clusters in the 
recovered Beers-Millennium catalogue. In the bottom panel 
of Fig. [TJ we show the scatter between the true and recov- 
ered halo masses in the Beers-Millennium catalogue. The 
scatter is large because of the considerable loss of informa- 
tion that occurs when going from the 3-dimensional to the 



1-dimensional velocity dispersion. This effect is especially 
severe for low mass clusters which have only a low number 
of detectable members. This means that that a large fraction 
of the massive Beers-Millennium Clusters, on which we fo- 
cus our analysis, have their masses and radii overestimated, 
i.e. are low mass clusters in the Millennium Catalogue. Con- 
sequently, the galaxy population in the cluster outskirts will 
be strongly contaminated by interlopers. 

We now examine the interloper fraction in our clusters 
in more detail. Any galaxy which is identified as a member of 
a Beers-Millennium cluster by the bi-weight algorithm, but 
does not actually reside within i?200,Mm of the corresponding 
true Millennium cluster, is defined as an interloper. In Fig. 
[2] we show the interloper fractions as a function of cluster- 
centric radius for clusters with M200 = 10 14 — 10 15 /i -1 Mq. 
The cluster center is defined to lie at the position of the 
central cluster galaxy. The interloper fraction in the cen- 
ter is very low, but rises steeply towards the outskirts, as 
expected. Since we have used a similar method as vdL07, 
our finding indicates that a similar contamination might be 
present in their stud^Q. Dashed lines show the fraction of in- 
terlopers which are central galaxies in the Millennium Cat- 
alogue. This fraction is very small, showing that most of the 
interlopers are in fact satellite galaxies; they either belong to 
the extended filaments surrounding the cluster (and are as 
such part of the fof cluster), or are part of infalling nearby 
groups. Note that as many as 40% of all satellite galax- 
ies in the SAM actually reside beyond 7?2oo,Mill of their re- 
spective central galaxy. These interlopers have by definition 
experienced environmental effects in the SAM and cannot 
be considered true "contamination" from the field; however, 
the time at which they fell into the cluster will differ from 
galaxies that reside within -R200 of the cluster center. 

We conclude that there is a high interloper fraction in 
clusters identified with a typical cluster-finding algorithm, 
and a relatively large scatter between true and recovered 
mass. This means that it is crucial to mimic the observa- 
tional cluster finding process when examining environmental 
effects in models. 



3 OBSERVATIONS 

3.1 Defining passive and active galaxies 

We use the specific star formation rate (SSFR) to distin- 
guish passive and active galaxies. The reason for our choice is 
that the SSFR is a physical quantity with a straightforward 
meaning. It can be directly taken from the SAM without 
any assumptions about stellar population synthesis models. 
However, the SSFR is not a directly observed quantity. 

Brinchmann et al. (2003) have derived SSFRs for galax- 
ies in the SDSS using emission lines, the D4000 spectral 
index and colours, with colours used for sampling star for- 
mation outside the region covered by the fiber. Salim et al. 
(2007) have redetermined SSFRs for a subset of these galax- 
ies from UV and optical photometry, and while they find 

1 It is however likely that the sample of vdL07 is biased towards 
higher cluster masses, since these are observationally easier to 
detect, which could make our contamination slightly higher than 
in vdL07. 
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Figure 1. Top panel: The fractional distribution of cluster masses 
for the 625 clusters in the vdL07 catalogue (black line) for the 
4393 clusters in the Beers-Millennium-Catalogue (blue dashed 
line) and for the 7586 clusters in the Millennium Catalogue (red 
dotted lines). Bottom panel: The relation between true and re- 
covered halo masses for the 4393 clusters in the Beers-Millennium 
Catalogue. Clusters with masses around 10 13 h —1 MQ in the Mil- 
lennium Catalogue suffer from a particularly large scatter in the 
halo mass, due to their low number of member galaxies which 
causes a large spread in ID velocity dispersion. 



good agreement for the "star-forming" class of the Brinch- 
mann et al. (2003) sample, they also find that star forma- 
tion rates are overestimated for the "star forming low S/N", 
the "Composite", the "AGN" class and the galaxies with- 
out measured Ha (see Salim et al. 2007 and Schiminovich 
et al. 2007 for more details). We therefore decide to use the 
SSFR from Brinchmann et al. (2003), but with a rough cor- 
rection, bringing it into statistical agreement with Salim et 
al. (2007). For "low S/N" galaxies, "AGN" and "composite" 
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Figure 2. Fraction of interlopers (i.e. galaxies that reside out- 
side -R200 in the Millennium Catalogue) in the Beers-Millennium 
clusters with masses 10 14 — 1O 15 /i -1 M0, in four different loga- 
rithmic stellar mass bins, as indicated. Note that the definition of 
an interloper used here is very strict and many galaxies classified 
as interlopers actually belong to the (non-spherical) large-scale 
structure associated with the cluster. Dashed lines show the frac- 
tion of interloper galaxies which are centrals in the Millennium 
Catalogue. 



galaxies, and galaxies with no Ha, we scale the Brinchmann 
et al. (2003) star formation rates down by 0.2, 0.4 and 1.0 
dex to correct for the offsets seen in Fig. 3 of Salim et al. 
(2007). The SSFR for the "star-forming" galaxies are used 
without any correction. 

In Fig. [3] top panel, we show the SSFR as a function of 
stellar mass for central galaxies in the SAM. A very clear star 
forming sequence at around log(SSFR)=-10 can be seen. In 
the bottom panel of the same figure, we show the corrected 
SSFR for central galaxies in the SDSS. These were deter- 
mined using the Yang et al. (2007) group catalogue and not 
vdL07, since it also includes low mass groups and thus pro- 
vides much better statistics. Every galaxy is weighted by 
1/Knax, where V ma x is the maximum volume out to which 
a galaxy of that magnitude would be detected in the sur- 
vey. We define "passive galaxies" as those with log(SSFR) 
< -11 both in the SAM and in the SDSS, because this cut 
corresponds roughly to the location of the minimum in the 
bimodal distribution of SSFR in both the model and the 
observations. 

Comparing the overall distributions in the SDSS and 
the SAM, we see that the shape of the star forming sequence 
is tilted in the SDSS, unlike in the SAM. This problem has 
been noted before (e.g. Somerville et al. 2008) , and seems to 
be generic to current SAMs (but see Neistein & Weinmann 
2009) 



3.2 Observational Results 

Here, we present the observational results that will be used 
as the basis for our model comparisons in what follows. 
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Figure 3. The SSFR as a function stellar mass in the SAM and 
the SDSS, for central galaxies only. Red lines indicate how we 
divide galaxies into "passive" and "active" in all of what follows. 
In the SDSS, we use the Brinchmann et al. (2003) SSFRs, cor- 
rected such that the offset with respect to Salim et al. (2007) is 
removed. In the SAM, galaxies with zero star formation rate have 
been randomly assigned SSFRs between -11.6 and -12.4, and the 
same has been done for SDSS galaxies which have measured SS- 
FRs, but with values < —12.4, since such low SSFRs tend to have 
large errors. In the SDSS, only galaxies with z < 0.1 are used, 
and a V max -correction has been applied. The number of galaxies 
has been normalized to the same value in both panels. 

These include the passive fraction of central galaxies as a 
function of stellar mass, the fraction of passive galaxies as 
a function of cluster-centric radius, and the distribution in 
SSFR for satellite and central galaxies. For all observations, 
we apply volume- weighting to correct for the apparent mag- 
nitude limit of the SDSS. 

In our Fig. we show as crosses with errorbars the 
passive fraction of central galaxies as a function of stellar 
mass defined using the group catalogue of Yang et al. (2007). 
As we show in Neistein & Weinmann (2009), similar results 
are obtained if the SSFRs of Salim et al. (2007) are used 
directly. 

We use the cluster catalogue of vdL07 to determine 



passive fractions as a function of projected cluster-centric 
radius in the SDSS (crosses with errorbars in Fig. [9}. All 
galaxies within 3 • o\n of the cluster center are used in this 
plot. Errorbars are determined using jackknife resampling. 
We have checked that results are virtually unchanged if we 
define volume-limited samples that are complete down to a 
given limiting stellar mass, as done by von der Linden et 
al. (2009). For the three lower stellar mass bins, the passive 
fractions as a function of cluster-centric radius are nearly 
identical to what has been found by von der Linden et al. 
(2009), which is somewhat surprising, since they have de- 
fined passive galaxies completely differently, with a method 
based on the light-weighted age of the stellar population, 
rather than the current SSFR as done here. This indicates 
that our way of classifying galaxies into passive and active 
is robust, at least at masses below log(M//i,- 1 M Q )=10.75. 
At higher masses, passive fractions are systematically lower 
by 10-15 % than in von der Linden et al. (2009). 

In Fig. [10] and Fig. [TT] we show as black histograms 
the distribution of SSFRs for SDSS central galaxies as de- 
fined by Yang et al. (2007), and for satellite galaxies de- 
fined here as galaxies located at a distance of less than -R200 
and within 3 • ctid from the center of a cluster with mass 
log(M/h~ 1 M Q ) > 14 in the catalogue of vdL07. A clear 
bimodality can be seen for both galaxy populations. 



4 DARK MATTER EVOLUTION OF 
SUBHALOES 

In this section, we investigate the stripping of dark matter 
subhaloes in the Millennium simulation. Stripping of sub- 
haloes has been studied numerically e.g. by Hayashi et al. 
(2003); Gao et al. (2004); De Lucia et al. (2004); Kazantzidis 
et al. (2004); Zentner et al. (2005); Giocoli et al. (2008); 
Diemand, Kuhlen & Madau (2007) and observationally us- 
ing weak gravitational lensing signals e.g. by Limousin et 
al. (2007) or Natarajan et al. (2009). The results presented 
here will be important for the new prescriptions for satellite 
evolution that will be introduced in the next section. 

In Fig. [4] we show the median fraction of dark matter 
which has been stripped since infall, for satellite galaxies in 
clusters and groups at z~0 (left hand and middle panel) and 
at z—1 (right hand panels), as a function of lookback time to 
infall. We have checked that results do not strongly depend 
on host halo mass down to 10 12 hT 1 Mq. Different linestyles 
denote different bins of Mj n f a n. Satellite galaxies are defined 
as galaxies which are part of a fof group, but are not the 
central galaxy of the most massive subhalo. For the satellite 
galaxies at z=0, we adopt two different methods for treating 
subhaloes which fall below the resolution limit (sometimes 
called "orphan galaxies"). In the left hand panel, we assume 
that their mass falls instantaneously to zero, while in the 
middle panels, we set set the mass equal to the effective 
resolution limit. For the highest mass subhaloes, the results 
are little affected by resolution and the curves in the left and 
middle panel are similar. For low mass subhalos, the two 
prescriptions (not surprisingly) give quite different results. 
The plots indicate that the Millennium simulation can be 
used to follow subhalo evolution to high accuracy for systems 
where Mi n f a n is larger than 10 11 h -1 Mq. 

Our results for such subhalos indicate that the dark 
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Figure 4. The median fraction of subhalo mass which has been stripped since infall as a function of the lookback time to infall. Only 
subhaloes containing galaxies with masses > 10 9 ' 5 /i — 1 Mq today in DLB07 are used. In the left hand panel, we show results for satellites 
at 2=0 if the mass of unresolved subhaloes is set to zero. In the middle panel, we show results for satellites at z=0 if the mass of 
unresolved subhaloes is set to the resolution limit of 1.7 X 10 10 h~ 1 Mq, for illustration. In the right hand panel, we show results for z=l, 
with the mass of unresolved subhaloes set to zero. 



matter is usually nearly completely stripped after the sub- 
halo has spent 5-7 Gyr in the cluster. Interestingly, we find 
that stripping is surprisingly similar for subhaloes spanning 
two orders of magnitude in dark matter mass at infall, if the 
dark matter mass of orphan galaxies is set to zero. We find 
that for satellite galaxies identified at z~0, roughly 15-20 % 
of the initial dark matter mass has been stripped per Gyr 
elapsed since infall. Note that the decline is not exponential. 
This means that stripping becomes more efficient the longer 
a subhalo has been a satellite. This can be explained by the 
fact that subhaloes which have already spent a significant 
amount of time as satellites have sunk to the center of the 
potential well (Gao et al. 2004), where stripping becomes 
fast, while a large fraction of subhaloes just having fallen 
in recently will be at the outskirts of the fof group where 
stripping is inefficient. Note that the stripping for an indi- 
vidual subhalo typically proceeds much less smoothly than 
the median shown here, as a large fraction of the stripping 
occurs at the first pericenter passage (e.g. Diemand et al. 
2007). 

For galaxies which are satellites at z=l (right hand 
panel of Fig. 3} , the dark matter stripping since infall has 
been significantly more efficient, with ~ 40 % of the initial 
dark matter mass stripped per Gyr in the median. This in- 
dicates that the stripping efficiency increases towards higher 
redshift, which is expected due to the decrease of dynami- 
cal times. This result is in qualitative agreement with the 
findings of Giocoli et al. (2008) and Tinker & Wetzel (2009). 
Note that both Natarajan, De Lucia & Springel (2007) and 
Maciejewski et al. (2009) have found indications that the 
subhalo masses found with SUBFIND are systematically low 
with respect to the true mass enhancement, or the masses 
found with more accurate 6D substructure finders. This ef- 
fect is most severe in the inner regions of clusters, where the 
background density of the parent cluster is high (Natarajan 
et al. 2007). This means that our estimates of dark matter 



stripping might be slightly too high, especially for subhaloes 
close to the cluster center. An alternative, and perhaps more 
accurate, method for estimating the amount of stripping in 
subhaloes would be to combine SAMs with analytical pre- 
scriptions based on higher resolution dark matter simula- 
tions, like presented in Taylor & Babul (2001). 



5 THE MODELS 

In this section, we present both the standard SAM of DLB07 
and our modifications of this model. 



5.1 The standard model 

Here we give an overview of how gas physics is modelled in 
central and satellite galaxies in DLB07 and repeat the basic 
equations (as given and explained in more detail in Croton 
et al. 2006 and references therein) for clarity. 

For central galaxies, gas surrounding galaxies comes in 
three phases: The cold gas, the hot gas and the ejected gas. 
Gas accreting onto the galaxy from the IGM is added to the 
hot phase. From this reservoir, it can cool down to the cold 
phase. Cooling is assumed to proceed differently in the "hot 
halo regime" (where r coo i < R v i T ) than in the "rapid cooling 
regime" (where r coo i > Rvii)- In the hot halo regime, cooling 
rates are given by 

'cool 

Vv IT 



"icooi = 0.5m h 



R 



2 

vir 



(4) 



7?vir is the virial radius of the dark matter halo, Vvix the 
virial velocity and r CO oi is the radius where the local cooling 
time t coo i is equal to the halo dynamical time. The local 
cooling time is given by 



^cool 



3 m p kT 
2p 9 (r)A(T, Z) 



(5) 
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with m p the mean particle mass, k the Boltzmann constant, 
and A(T, Z) the cooling function according to Sutherland & 
Dopita (1993) as a function of the temperature of the gas, 
which is assumed to be equal to the virial temperature of 
the halo, and the gas metallicity. p g (r) is the hot gas density, 
which is assumed to have an isothermal profile 



Pg{r) 



47T_Rvirr- 2 



(6) 



If fcooi > Rvir, the cooling rate is set roughly to the rate at 
which new diffuse gas is added to the halo (see discussion in 
Croton et al. 2006). 

The energy input due to radio mode AGN feedback ac- 
cording to Croton et al. (2006) is assumed to partially offset 
the cooling in high mass haloes, giving rise to a modified 
cooling rate 



m* CO ol,mod m, coo l 



Lbh 

2~ 

vir 



0.5K 



with the mechanical heating energy generated by the black 
hole accretion Lbh given by 



Lbh = T]rn,BB\c 



(8) 



rj =0.1 is the standard efficiency with which mass is assumed 
to produce energy near the event horizon, and c the speed 
of light. The quiescent growth rate of the black hole is given 
by 

Kir N 3 



(9) 



, "IBH , , /hot \ , 

mBH = KAGN( 10^o )( 0T )( 200kms- 1 

where dibh is the black hole mass, /hot the fraction of the 
total halo mass in the form of hot gas and kagn a free 
parameter set to 7.5 • 10~ 6 Moyr _1 in DLB07. 
The cold gas then forms stars according to 

m-star = "SF (m co ld — m cr lt )/tdyn,disk (10) 

with £dyn,disk = fdisk/Kir, r di s k = 3( \/\f2)R v - 1T (motivated 
by the studies of Mo et al. 1998 and van den Bergh et al. 
2000), and A the spin parameter of the dark matter halo in 
which the galaxy resides. a$F is a tunable parameter and set 
to 0.03 in DLB07. The calculation of the critical gas mass 
follows Kauffmann (1996) and is based on the observations 
of Kennicutt (1998) of a threshold gas density below which 
stars do not form anymore: 

Kit 



merit = 3.8- 10 a ( ; 



-)( 



»"disk 



)M G 



(11) 



^OOkms- 1 "lOkpc' 

Supernova feedback heats cold gas back to the hot phase 
with an efficiency directly proportional to the star formation 
rate: 



e d isk ■ Am, 



(12) 



based on the observations by Martin (1999), with edisk = 
3.5 set by observational data. Gas in the hot phase can be 
transported to the ejected phase due to the excess energy 
present in the hot halo after reheating, with an efficiency 
inversely proportional to the depth of the dark matter halo 
potential: 



Am,ejected = (dialo 



V 2 



e d isk)Am» 



(13) 



with Vsn=630 km s 1 the mean energy in supernova ejecta 
per unit mass of stars formed, based on a standard IMF and 



standard supernova theory, and ehaio=0.35 tuned to repro- 
duce observations. Per halo dynamical time, half of the gas 
in the ejected phase is assumed to be reaccreted to the hot 
phase: 



A?71 C j Cc tcd — "/ ' ^T^cjected/^dyrj 



(14) 



with 7=0.5 and t^yn the halo dynamical time. 

In Fig. [5] we plot the fraction of the total gas in the 
cold, hot and ejected phase in central galaxies at 2=0 both 
as a function of galaxy stellar mass and of dark matter halo 
mass in the DLB07 SAM. Clearly, the ejected phase strongly 
dominates both at low stellar and halo masses, where the hot 
phase is relatively unimportant. This is a consequence of the 
rapid cooling and the very efficient SN feedback in low mass 
haloes, and indicates that the ejected phase, and how it is 
reaccreted, is of crucial importance for the evolution of low 
mass galaxies (see also Oppenheimer et al. 2009). 

Satellite galaxies in DLB07 are defined as all galaxies 
which are member of fof groups, but not the central galaxy 
of the most massive subhalo. In all our models presented 
below, we will follow this definition. If a galaxy becomes a 
satellite in DLB07, both the hot and ejected gas is removed. 
Any hot and ejected gas produced after infall is stripped 
as well. Due to the strong efficiency of SN feedback in the 
model, this leads to a quick depletion of cold gas and to 
the cessation of star formation. All stripped satellite gas is 
added to the hot gas of the central galaxy. 



5.2 Changes to the standard model 

In what follows, we describe our basic changes to the stan- 
dard DLB07, which are ingredients of all our models 1) - 4) 
described below. 

• We allow cooling from the hot to the cold gas, and rein- 
corporation from the ejected to the hot gas, for satellite 
galaxies. More details are given below. 

• In the DLB07 SAM, processes like cooling and star for- 
mation in galaxies are correlated with the properties of the 
corresponding halo, as apparent from the equations given 
above. After a galaxy becomes a satellite, this is not neces- 
sarily appropriate. For example, it might not generally be 
true that the disk radius of a satellite decreases in propor- 
tion to the radius of its dark matter halo, as assumed in 
DLB07. De Lucia & Helmi (2008) therefore fix the disk ra- 
dius of the satellite at infall. It is also unclear how the hot 
gas halo redistributes itself after part of the dark matter has 
been stripped. In our modified model, we fix R v i r , A and V v i T 
for satellites in eq. l4H4l at infall and do not let it evolve any- 
more. This means that star formation, cooling, ejection and 
reincorporation for satellites are calculated according to the 
subhalo properties at infall. We still follow the true evolu- 
tion of the dark matter subhalo in detail, as it is used for (i) 
determining dynamical friction timescales and (ii) our new 
recipes for satellite stripping, as described below. 

• Unlike DLB07, we allow AGN feedback for satellite 
galaxies as described in eq. [7] [5] and [5] in our models 1) - 
4), again using R v i T and V V - 1T as determined at infall. As ra- 
dio AGN activity has been observed in satellite galaxies (e.g. 
Best et al. 2007; Pasquali et al. 2009a), this step seems well 
justified. 

• It has been found in previous work (e.g. Weinmann et 
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Figure 5. The average fraction of gas in the cold, hot and ejected 
phase, at 2=0, as a function of stellar mass (top panel) and of 
halo mass (bottom panel), for central galaxies with stellar masses 
above log(M//i- 1 M Q )=9.5 in the SAM of DLB07. 



al. 2006b) that satellite galaxies in DLB07 are too red. As 
mentioned before, DLB07 assumes that all hot and ejected 
gas is removed from satellites. We find that letting satellite 
galaxies keep their entire hot gas, while completely stripping 
the ejected phase, still produces too many passive low mass 
satellites, as only a small fraction of the gas in these galaxies 
is actually in the hot mode (see Fig. [5} . We therefore use a 
different approach. At each timestep, we assume that some 
fraction / str i P of Af diffusc = Af hot + M cjcctcd is stripped. We 
assume that the ejected gas is always stripped first, which 
means that the hot gas reservoir is only stripped when /strip • 
Afdiffusc > Afojoctcd- All the stripped gas is always added 
to the hot gas of the central galaxy in the corresponding 
FOF group. In our models 1) - 3), we strip the diffuse gas 
halo of satellite galaxies in proportion to the dark matter 



subhalo. This means that if a given satellite galaxy loses 
a fraction of its dark matter mass between two subsequent 
timesteps i-1 and i, the same fraction /strip, i of the current 
mass in diffuse gas is removed. 

Mdma 



f, 



strip.; 



= l 



M E 



(15) 



To calculate cooling rates, we need to make an assumption 
about the hot gas density profile of satellites. Here we sim- 
ply assume that the the remaining hot gas redistributes it- 
self according to eq. [5] (with R v i T as recorded at infall) after 
each stripping event. This means that the gas becomes pro- 
gressively more diluted, which additionally decreases cool- 
ing rates. This method differs from Font et al. (2008), who 
determine a stripping radius and remove all gas outside of 
this radius, causing the hot gas halo to become increasingly 
compact. We think that our simpler approach is justified, 
because the distribution of the hot gas in satellite galaxies 
is highly uncertain; continued energy input by SN feedback 
might well cause the hot gas halo to expand after a strip- 
ping event has occurred. Also, is is not entirely clear how 
-Rvir and V v i t should be calculated for subhaloes. 

Finally, we note that if the satellite galaxy falls below the 
resolution limit, all the hot and ejected gas is stripped at 
the timestep of transition. After that, the gas content of the 
galaxy only evolves as a result of star formation and SN 
feedback. 



5.3 Specific Modifications 

We carry out a series of progressive modifications to the 
semi-analytic models of DLB07, denoted modifications 0) - 
3), as described below. We also try out a simpler approach, 
which does not follow the stripping of the DM haloes in 
detail, denoted modification 4). 

• In our modification 0), the only change with respect to 
DLB07 is that we use _R v i r and Kdr as recorded at infall in 
eq. 141 14l as described in section [531 second point. Since all 
hot and ejected gas is still removed from satellites, cooling 
and reincorporation rates for satellites are still zero, as in 
DLB07. We have checked that results are virtually indistin- 
guishable from DLB07, with the only exception of satellite 
morphologies, which are brought into better agreement with 
observations. We plan to study this issue in future work. 

• In our modification 1), we implement all the changes 
described in section 1531 i.e. we strip the diffuse gas in pro- 
portion to the dark matter. 

• For modification 2), we additionally decrease the qui- 
escent hot gas accretion rate to the black hole kagn (as 
defined in eq. [§]) by a factor of 5, and we let satellite galax- 
ies residing in host haloes with masses below 10 12 /i _1 M© 
keep all their hot and ejected gas. The motivation for this is 
discussed in section 16.31 For illustration, we also present a 
modification 2b), which only includes the decrease in kagn- 

• While the dark matter subhalo mass for most simulated 
satellite galaxies continuously decreases, there are also cases 
in which a subhalo seems to experience dark matter growth. 
The reason for this is mainly that the SUBFIND algorithm 
(Springel et al. 2001), which is used to identify subhaloes, 
defines the boundary of a subhalo as the point where the 
density is equal to the background density. The same sub- 
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halo will therefore have a lower mass assigned if it is close to 
the center of its host halo, compared to when it is residing 
in the outskirts. This can lead to artificial mass fluctuations 
and thus to overstripping of diffuse gas in our approach. We 
correct for this effect in our modification 3). If the dark mat- 
ter subhalo mass at a given timestep is higher than in the 
previous timestep, we add the appropriate amount of gas to 
the hot gas halo of the satellite. 

• Modification 4) differs from the previous ones, since it 
does not rely on following the stripping of dark matter sub- 
haloes. This prescription may be used by SAMs which do 
not follow the stripping and the orbits of subhaloes after 
they have become part of a fof group (like e.g. Somerville 
et al. 2008). We simply store the mass in diffuse gas at the 
time of infall, and strip some fixed fraction of this initial 
(not the current!) diffuse gas mass per Gyr. We determine 
this fraction by tuning the model to fit observations. We 
also require that the median stripped dark matter fraction 
roughly agrees with the results calculated directly from the 
simulations, which have been discussed in section [4] We find 
that a good solution is to strip 20 % of the initial gas mass 
per Gyr if the subhalo mass at infall is < 5 x 10 u ft _1 M Q) 
and 10 % per Gyr if it is above this mass. As in our models 
2) and 3), we decrease the AGN feedback efficiency, satel- 
lites residing in groups with masses below 10 12 /i _1 M© arc 
not stripped at all, and as in models 1) - 3), we always first 
strip the ejected, and then the hot gas. 

A summary of our models is given in Table 1. 



6 RESULTS 

In this section, we compare the model results to the obser- 
vational results described in section I3T21 For all plots in this 
section, we use the same subvolume of the Millennium Sim- 
ulation as used in section [2~4l The bi- weight algorithm de- 
scribed in the same section is re-applied for each new model. 
For the model results shown in the plots in this section, cen- 
tral galaxies are defined as the galaxy belonging to the most 
massive subhalo in a given FOF group, while satellites are de- 
fined to be those galaxies located at a projected distance of 
less than -R200 and within 3 • am from the center of a cluster 
with log(M//i _1 M ) > 14 (with M determined using the 
bi- weight algorithm). 

6.1 Modification 

Implementing modification 0) gives results that are virtu- 
ally indistinguishable from DLB07. In Fig. [SI we compare 
the passive fraction of central galaxies in modification 0) 
(blue dotted line) with observations (data points). Surpris- 
ingly, the agreement is rather poor. At intermediate stellar 
masses, the passive fraction is underproduced, while it is 
overproduced at high stellar masses. The disagreement at 
intermediate stellar masses might partially be due to the 
contamination of the central galaxy sample of Yang et al. 
(2007) by satellite galaxiefl. However, it is more likely that 

2 The fraction of galaxies which are satellites at masses 
log(M st ar//i _1 M©) ~ 10 in is higher by ~ 15 % in the SAM 
than in Yang et al. (2007). Under the rather extreme assump- 



Table 1. Overview of the five modified version of DLB07 dis- 
cussed in this paper. Modifications 1) - 4) include all the changes 
discussed in section r5.2l modification 0) only the first point, kagn 
is the the quiescent hot gas accretion rate to the black hole in 
Mq/yv, as defined in Croton et al. (2006). For modification 4), 
/1 and /2 are the fractions of the initial diffuse gas stripped per 
Gyr. fi is for an initial subhalo mass of < 5 ■ lO 11 /l- 1 M , h for 
higher masses. 
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DLB07 underproduces the number of passive central galax- 
ies at intermediate masses. This problem has already been 
noted by Baldry et al. (2006) and Bolzanella et al. (2009) 
who found that Croton et al. (2006) model did not produce 
enough blue central/field galaxies. The apparent overpro- 
duction of passive galaxies with low SSFR at high stellar 
masses, on the other hand, is more puzzling, since the frac- 
tion of red, massive galaxies seems to be in agreement with 
observations (e.g. Weinmann et al. 2006b). There are two 
potential explanations for this apparent discrepancy. First, 
precise determination the SSFRs of massive, red galaxies is 
difficult, and it is entirely possible that they may have been 
overestimated in our data. We note that von der Linden et al. 
(2009) find higher passive fractions at \og(AI staT / h" 1 Mq) ~ 
10.7 than we do here. Second, it is important to remember 
that passive fractions defined according to colour are dif- 
ferent from those defined according to SSFRs. Colours are 
sensitive to dust attenuation as well as the fraction of young 
stars in the galaxy, so the agreement with the fraction of 
red objects found in previous work (e.g. Croton et al. 2006) 
may just have been fortuitous. The main goal of the current 
paper is to improve the treatment of environmental effects, 
so we will not investigate this issue in detail; we simply de- 
crease AGN feedback in some of our models, as described 
below. 

In Fig. [9] we compare modification 0) (blue dotted line) 
with the passive fractions as a function of cluster-centric 
radius. Clearly, the model overproduces passive satellites, 
in agreement with previous findings (e.g. Weinmann et al. 
2006b). In Fig. [TOl and Fig. [11] we show the distribution 
of SSFR for central and satellite galaxies in modification 
0) (blue lines) compared to observations (black lines). The 
number of satellite galaxies with high star formation rates is 
clearly too low compared to the observations. For the cen- 
tral galaxies, the location of the active peak is not correctly 
reproduced, as a result of the missing tilt in the relation be- 
tween SSFR and stellar mass discussed before. For all models 
discussed below, the location of the two peaks for the central 



tions that the satellite fraction in the SAM is correct, that these 
additional 15 % of satellites are wrongly counted as centrals in 
the Yang et al. (2007) group catalogue and that all of them are 
passive, we would be able to roughly account for the observed 
difference. 
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galaxies is very similar, and only the relative heights change. 
In what follows, we will therefore simply concentrate on the 
passive fractions for the central galaxies. 

6.2 Modification 1 

In our modification I), we strip the diffuse gas at the same 
rate as the dark matter of subhaloes is stripped due to tidal 
effects. We find that the overproduction of passive central 
galaxies get slightly worse (see Fig. solid cyan lines). The 
reason for this is that less hot gas is now available for cooling 
to the central galaxies. For the satellite galaxies, results are 
improved with respect to DLB07 in the two lower mass bins, 
while results are similar in the two higher mass bins (see 
Fig. [9] solid cyan lines). In Fig. 1111 cyan lines indicate the 
distribution of SSFRs for satellite galaxies. 

For illustration, we compare the star formation, dark 
matter and gas accretion histories of two random satellite 
galaxies at z—Q with \og{M st a.T: / h" 1 Mq) ~ 10. Results are 
shown for the DLB07 prescriptions (top panels) and for 
modification 1) (bottom panels) in Fig. [6] In the left hand 
(right hand) panels, the galaxy becomes a satellite at around 
a lookback time of 3 Gyr (9 Gyr) , which causes the hot and 
ejected gas reservoir to fall to zero in DLB07. In our mod- 
ification 1), stripping is much slower. Despite the fact that 
we always strip the ejected gas first, the hot gas reservoir 
is depleted more quickly for the galaxy in the right hand 
panel. This is due to the fact that nearly all hot gas that is 
added by SN feedback or that is re-incorporated from the 
ejected phase cools down efficiently in low mass galaxies. 

In Fig. we show the relation between the median 
SSFR and the fraction of the dark matter subhalo mass 
stripped since infall (top panel), and the lookback time to 
infall (bottom panel) in modification 1) for galaxies with 
stellar masses log(M / 7i _1 Mq) = 10 — 10.5. It can be seen 
that galaxies typically only become passive after around 80 
% of the dark matter subhalo (and thus of the diffuse gas) 
has been stripped, and only after around 5 Gyr since infall. 

6.3 Modification 2 and 2b 

We try to address the overproduction of passive central 
galaxies caused by modification 1) by lowering kagn, as de- 
fined in eq.(9] Green dot-long dashed lines in Fig. [8]show the 
effect of modification 2b), in which kagn is lowered by a fac- 
tor of 5. The satellite stripping is treated exactly as in model 

1) . While results are improved at the massive end, the pas- 
sive fractions fall well below the observations at intermedi- 
ate stellar masses. After some experimenting, we found that 
this problem can largely be solved if satellite galaxies in host 
haloes with masses below W^h -1 Mq keep all their hot and 
ejected gas. The reason is that this substantially reduces the 
amount of gas available for cooling to the central galaxy in 
such haloes. We implement this change in our modification 

2) , shown as green dashed lines in Fig. [8] Agreement with 
observations is improved compared to DLB07, although it 
is still not perfect. 

Why should gas stripping cease to be effective in halos 
with masses below 10 12 Mq? One reason is that a hot gas 
halos no longer form at these halo masses (see Fig. , which 
means that ram-pressure stripping will no longer be effec- 
tive. Our results are in agreement with McGee et al. (2009), 
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Figure 7. We show the median SSFR and the region where 68 % 
of the values lie, as a function both of the fraction of mass which 
has been stripped since infall (top panel) and as a function of 
the lookback time to the last time the galaxy became a satellite 
(bottom panel). Results here are shown for galaxies with masses 
of log(M//i- 1 M )=lO-lO.5, at z=0, in our modification 1). 

who find that observations are best reproduced under the 
assumption that no significant quenching of star formation 
occurs in host haloes with masses below ~ 1O 12 /i _1 M . We 
note that significant dark matter stripping does still occur 
in halos with masses below 1O 12 /i _1 M . Our results should 
thus be taken as an indication that the prescription for gas 
stripping that we have introduced, in which the gas is re- 
moved at the same rate as the surrounding halo loses mass, 
is simply an approximation that appears to give roughly the 
right answer for satellite galaxies in massive halos. It should 
not be regarded as a correct representation of the physics 
that controls the stripping of the gas. 

Modification 2) is quite successful in reproducing pas- 
sive fractions as a function of cluster-centric radius (see Fig. 
green long-dashed lines). The remaining disgreement is 
caused by the fact that massive central galaxies are still 
slightly too passive in this model - this problem propagates 
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lookback time [Gyr] lookback time [Gyr] 

Figure 6. Here we show the histories of two random galaxies with ]og(M s ta.v/h~ Mq) ~ 10, for the DLB07 model (top panels) and 
in our modification 1) (bottom panels). We plot the masses in cold, hot and ejected gas and of the DM subhalo, and log(SSFR) (with 
values of the latter given at the right axis of the right hand plots). For both example galaxies, the hot and ejected gas mass drops to 
zero after infall in DLB07, which happens around 3 Gyr before present in the galaxy shown in the left hand panels, and at around 9 
Gyr before present in the galaxy shown in the right hand panels. The galaxy in the right hand panel shortly becomes a central again at 
around 7 Gyr. At 2=0, the galaxy in the left hand panel has a stellar mass of log(M//i- 1 M Q )=10.1 (10.07) in our modification 1) (in 
DLB07). The galaxy on the right hand panel has a mass of log(M/fc. _1 M ) 10.18 (10.09). 



to the satellite galaxies irrespective of the stripping that we 
implement. However, if we look at the detailed distribution 
of the SSFR of satellite galaxies (Fig. 1111 green histogram), 
we find that the star forming peak of the satellite galaxies is 
somewhat too low, and that that there are too many satellite 
galaxies with intermediate SSFR present. 



6.4 Modification 3 

In modification 3), we correct for fluctuations in subhalo 
masses, by letting satellite galaxies "accrete" gas from the 
hot gas reservoir of their host halos, if their dark matter 
subhalo experiences apparent mass growth. While we find 
that the distribution of SSFR of satellites is now better 
reproduced (Fig. 111! red histogram), the passive fractions 
(Fig. [9] red dashed lines) do not agree as well with the data. 
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Figure 8. Here we show the impact of our different modifications 
to the SAM on the passive fraction of central galaxies. Different 
lincstyles correspond to different modifications, as indicated and 
described in the text. Crosses with errorbars show the observa- 
tional results. 



6.5 Modification 4 

By construction, modification 4) produces results that are 
in good agreement with observations, both in terms of the 
distribution of SSFR and the passive fractions as a function 
of cluster-centric radius. The results are shown in Fig. [8] - 
1111 as magenta dot-dashed lines. Note that it is not possi- 
ble to obtain a similar level of agreement if we assume that 
gas is stripped exponentially. While it is relatively easy to 
reproduce passive fractions, we find that an exponential de- 
cline of the gas reservoir cannot at the same time produce 
a bimodal SSFR distribution for the satellites. Also Balogh 
et al. (2009) have pointed out that reproducing the detailed 
bimodal colour distribution of satellite galaxies with SAMs 
is not trivial. 



6.6 The impact on the global properties of the 
galaxy population 

In this section, we investigate the properties of the global 
galaxy population in our alternative models. We have not 
tuned our models to match these and it would thus not be 
surprising if the models were not to agree with the obser- 
vations. Here, the full simulation box is used. Note that for 
all results shown here, there is no visible difference between 
DLB07 and our model 0). 

The stellar mass functions in the different models are 
compared to the observations of Li & White (2009), cor- 
rected according to Guo et al. (2009), in the top panel of 
Fig. 1121 All models overproduce the massive end of the stel- 
lar mass function. In our models 2) - 4), this problem is 
accentuated. This is due to the decreased efficiency of AGN 
feedback that we need in order to reproduce the fraction of 
passive central galaxies. This clearly indicates that simply 
decreasing AGN feedback is not a viable solution; maybe, 
it would be more appropriate to make AGN feedback more 
stochastic, for example by shutting down gas cooling after 




0.05 



-12 -11 
log(SSFR) [yr 



Figure 10. The distribution of the central SSFR in DLB07 com- 
pared to the SDSS central galaxies according to Yang et al. (2007). 
Galaxies with zero SSFR in the SAM are randomly distributed 
between log(SSFR)=-11.6 and log(SSFR)=-12.4, and the same is 
done for SDSS galaxies with measured log(SSFR) < -12.4. Re- 
sults for our alternative SAMs are not shown since the location 
of the peak is identical with DLB07, and only the relative peak 
heights change, which is information already given in Fig. [8] 



major mergers (e.g. Hopkins et al. 2008). Interestingly, the 
overproduction of luminous galaxies is more severe in our 
modification 3) than in our modification 2), which only dif- 
fer in details of the implementation of satellite gas stripping. 
This indicates that the treatment of satellite can have a sig- 
nificant impact on the general galaxy population. 

In Fig. 1131 we compare the total bj-band luminosity 
function (LF hereafter, top panel) and the red 6j-band LF 
(bottom panel) for the different models to observations. 
Galaxies are classified as red if B — V < 0.8, following Cro- 
ton et al. (2006). Observations are taken from Norberg et al. 
(2002) and Madgwick et al. (2002). Again, the impact of de- 
creasing AGN feedback in our modifications 2) - 4) is clearly 
apparent. In the bottom panel of Fig. 1131 it can be seen that 
all models shown here have some problems with reproducing 
the red LF. Interestingly, our improved treatment of satellite 
stripping does not entirely solve the overproduction of faint 
red galaxies noted by Croton et al. (2006). We note that our 
ability to link gas stripping to the evolution of the dark mat- 
ter mass of the subhalo will break down near the resolution 
limit of the simulation. Right now we assume that all the 
gas is stripped when the subhalo is finally disrupted. This 
will likely result in effective gas stripping efficiencies that 
are too high for the low mass satellites. It thus might be ad- 
visable to use a different method for galaxies that reside in 
subhalos with masses at infall close to the resolution limit. 
One possibility is to adopt modification 4) in this regime. In- 
deed, Figure 11 shows that modification 4) gives improved 
agreement with the faint end of the red LF. 

Another problem becomes apparent at intermediate lu- 
minosities in our modifications 2) - 4), where the number of 
red galaxies is too low. This may be related to the decrease 
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Figure 9. Here wc show the impact on our modifications on the passive fraction as a function of cluster-centric radius in four different 
stellar mass bins, as indicated on the top of the panel. Different linestyles correspond to four alternative SAMs, as indicated and described 
in the text. Crosses with errorbars show the observational results. Only clusters with masses above lO 14 h -1 Af0 are used. 



in AGN feedback. Stronger dust attenuation in these galax- 
ies may solve this problem, but we do not investigate this 
issue in more detail here. 

Our results show that the implementation of the strip- 
ping of the satellite galaxy gas supply in proportion to the 
dark matter does not impact the global properties of galax- 
ies very significantly, except by decreasing the abundances of 
faint red galaxies to some degree. However, decreasing the 
AGN feedback, and shutting off satellite stripping at halo 
masses below 10 12 ft -1 Mq in order to better reproduce the 
passive fraction of central galaxies clearly leads to inconsis- 



tencies with the LF. It is likely that in order to simultane- 
ously reproduce the passive fraction of centrals, the LF of 
red galaxies and the bright end of the LF, a more careful 
retuning of the model, including the prescriptions for AGN 
feedback and dust attenuation, will have to be carried out. 
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Figure 11. The distribution of the satellite SSFR in DLB07 and our 4 modifications compared to observed satellite SSFRs. Galaxies 
with zero SSFR in the SAM are randomly distributed between log(SSFR)=-11.6 and log(SSFR)=-12.4, and the same is done for SDSS 
galaxies with measured log(SSFR) < -12.4. Both in the models and the observations, satellite galaxies are defined as galaxies residing 
within -R200 from the centers of clusters with masses > 1O 14 /i _1 M0 as found by the bi-weight algorithm described in section [2.41 



7 DISCUSSION 

7.1 Treatment of satellite galaxies in the SAM 

Recently, several observational and theoretical studies have 
shown that the removal of the gas reservoir in satellite galax- 
ies is overly efficient in most current SAMs (McCarthy et al. 
2008; Simha et al. 2008; Font et al. 2008; Weinmann et al. 



2006b; Wang et al. 2007). In this work, we present a simple 
way to treat environmental effects. 

It seems physically plausible that the stripping effi- 
ciency may depend on the mass of the host halo, on the 
mass of the satellite galaxy (e.g. Bekki 2009), on orbital 
parameters (e.g. McCarthy et al. 2008) and on the age of 
the Universe (e.g. Giocoli et al. 2008) These dependencies, 
however, are not well constrained. It is not clear whether 
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Figure 12. The stellar mass functions in the standard model of 
DLB07 (blue dotted lines) and in our modifications of this model, 
as indicated. The observed mass function is a modified version of 
the one presented in Li & White (2009), corrected according to 
Guo et al. (2009). 



hydrodynamical simulations can properly predict the distri- 
bution of the (probably multiphase) hot and ejected gas, or 
deal with the presence of shielding magnetic fields. Smooth 
particle hydrodynamics codes may have difficulties dealing 
with dynamical instabilities occurring in regions with steep 
gas density gradients (e.g. Agertz et al. 2007). 

In this work, we suggest using the approximation that 
the diffuse gas halo of the satellite (which includes the hot 
and ejected gas) is stripped at the same rate as the dark 
matter halo. This is a very simple model which does not 
treat the detailed physics of gas stripping, but is easy to 
implement. We find that observations at z~0 can also be 
well reproduced by a simple recipe in which 10-20 % of the 
initial gas reservoir in the halo is stripped per Gyr. We stress 
that in order to reproduce observations, the the diffuse gas 
must be stripped linearly, and not exponentially. 

It is not yet clear whether environmental effects at 
higher redshift are correctly reproduced by this method. We 
find that the dark matter stripping seems to be more efficient 
at high redshift than at z~0. It will be a natural next step 
to compare our models to observations of galaxy properties 
as a function of environment at higher redshifts. 

Eventually, it will be important to understand the na- 
ture of the gas that is found around satellite galaxies. "Star- 
vation" is usually defined as the stripping of the hot gas 
halo (Larson, Tinsley & Caldwell 1980). However, below a 
mass limit of around W 12 h -1 Mq, gas cools rapidly and thus 
does not stay in a hot gas reservoir for a significant amount 
of time. Consequently, most gas in low mass haloes resides 
in the so-called "ejected" , and not in the "hot phase" , with 
the ejected mode composed of gas that has been expelled 
from the galaxy by supernovae. The state of this ejected gas 
is very poorly understood: it could consist of cold clumps, 
or it might form an extended disk, or it might have high 
temperatures and low densities. The only certainty we have 
about this phase is that it must remain bound to the galaxy 



„ -2 - 



CD 

6 

CO 

I 

O 

Oh 

S 

CO 



-4 - 



-6 - 



-8 



**xx 1 1 1 1 


— 


— 






w - 








X ', 




X 




x* - 


- x Norberg et al. 




. 1 i i . . 1 


1 1 1 



-15 



-20 



M b - 51ogh 
red LF 



00 
CCS 

6 

CO 

o 
a 

CO 



-4 - 



-6 



-8 



t — r 



T 



~\ — r 



x x 



mod. 1) 
mod. 2) 
mod. 3) 
mod. 4) 
DLB07 

Madgwick et al. 02 




J I I L 



1 



J L 



-15 



-20 



M b - 51ogh 

Figure 13. The total 6j-band luminosity function (top panel) 
and the fej-luminosity function of red galaxies alone (bottom 
panel) for the standard model of DLB07 (blue dotted lines) and 
our modifications of this model, as indicated. Observations are 
taken from Norberg et al. (2002) and Madgwick et al. (2002) and 
are shown as crosses and without errorbars. The split into blue 
and red galaxies has been done at B — V = 0.8. 



and form a significant part of the reservoir for future star 
formation in low mass systems. At low stellar masses, strip- 
ping of this ejected phase is the main driver of environmental 
effects in the local Universe. 

In this work, we have not considered the abundances of 
satellite galaxies as a function of halo mass or the correlation 
function on small scales; in order to reproduce these observa- 
tions, it is likely that tidal disruption and satellite-satellite 
mergers have to be taken into account, as SAMs tend to 
overproduce the abundance of satellite galaxies (e.g. Liu et 
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al. 2009). Tidal disruption might also help with reproducing 
the bright end of the luminosity function, as the number of 
mergers would be reduced, and would explain the presence 
of significant intracluster light (Zibetti et al. 2005). Sugges- 
tions of how tidal stripping might be implemented into SAM 
can be found in Kim et al. (2009), Wetzel & White (2009) 
and Henriques & Thomas (2009). Also, we have not consid- 
ered the stellar ages and metallicities of satellites, for which 
Pasquali et al. (2009b) have found some interesting discrep- 
ancies between SAMs and observations. These may help to 
further constrain the physics of satellite stripping in future 
work. 



7.2 Comparison to a simple model of 
ram-pressure stripping 

Gunn & Gott (1972) were the first to suggest that the ram 
pressure of the hot intracluster medium may strip gas from 
galaxies. Ram-pressure stripping of the cold disk gas has 
since received considerable theoretical attention (e.g. Abadi 
et al. 1999; Schulz & Struck 2001; Roediger & Briiggen 
2007). Ram-pressure stripping of the hot gas halo, how- 
ever, has been investigated by only a few studies (e.g. Mori 
& Burkert 2000; Hester 2006; McCarthy et al. 2008). Mc- 
Carthy et al. (2008) find that ram-pressure stripping is more 
efficient than tidal stripping in depleting the hot gas halo of 
most satellite galaxies. Our simple prescription described in 
section [5] does not explicitely model this effect. In this sec- 
tion, we test whether we can reproduce observations by im- 
plementing ram-pressure stripping in a physically motivated 
fashion. 

Our model is based on the model of ram-pressure strip- 
ping by Font et al. (2008), but with a few important modi- 
fications. We follow the orbits of the satellites in detail, and 
allow stripping at every timestep, while Font et al. (2008) 
only allow stripping at infall and at each time the host halo 
doubles in mass. We assume that the ram-pressure stripping 
of the hot gas proceeds from the outside in, and does not 
affect the central density and hence the cooling efficiency of 
the satellite hot ga^f]. We allow the assumption of outside-in 
stripping to break down if the truncation radius of the hot 
halo reaches some fraction /? of the virial radius of the dark 
matter halo that hosts the satellite at the time of infall. At 
this point, we allow the density profile of the hot gas to relax 
to an isothermal distribution out to the virial radius -Rvir, 
and we strip all the ejected gas. 

We describe the detailed implementation of the model 
in the L-galaxies code. While -Rtrunc > /3R v i T , we describe 
the density profile of the hot gas for satellites with a modified 
version of eq. [5] following Font et al. (2008) : 



Pg(r) 



m hot+ m strip 
47rfl vir r 2 





if r < Rtrun 
if r > Rtrun 



(16) 



where m str ip is the sum of all hot gas which has been stripped 
from the galaxy since it became a satellite. -Rtrunc is set to 
i?vir when a galaxy first changes from central to satellite. 



3 Unlike Font et al. (2008), we do not limit the cooling radius 
(^cool, see eq. |4] and O by -Rtrunc In eq. [4] we replace m no t by 
m hot + m strip for satellite galaxies. 



(Note that 7? v ir and V v it for satellites are fixed at infall and 
then remain constant.) 

The stripping radius i? s tri P is calculated according to 
Gunn & Gott (1972) and McCarthy et al. (2008) as the 
radius where the ram pressure exceeds the pressure induced 
by the self-gravity of the satellite. For a singular isothermal 
sphere of hot gas, we determine R s tri P with 



a/5g,sat (^strip)Kdr,s 



cn(R)V 2 



(17) 



with R and V indicating the distance of the satellite to the 
cluster center, and its velocity relative to it. For the pressure 
exerted by the self-gravity of the satellite, we use a prefactor 
q = 2, following McCarthy et al. (2008). However, we note 
that removing this prefactor has very little impact on the 
results. If i? s trip < i?tranc, we remove all the mass enclosed 
between those two radii, and set -Rtrunc = .Rstrip. 

If -Rtrunc < /3-Rvir, We reset m st rip = and -Rtrunc = Rvir, 

which means that the hot gas gets diluted out to -Rvir, lead- 
ing to strongly decreased cooling, and increased vulnerabil- 
ity to stripping. Galaxies often fluctuate between the central 
and satellite status. Whenever a satellite galaxy becomes 
a central again, its R v i r is updated according to its cur- 
rent dark matter subhalo. Rtrunc is set to R v ir- After each 
stripping event, the hot gas mass in the satellite will be 
modified by cooling, ejection, SN feedback and reincorpora- 
tion. Before we calculate the stripping occurring in the next 
timestep, we reset the truncation radius such that 



M(r < -Rtrunc) = M h 



(18) 



Up to now, we have only discussed the ram-pressure 
stripping of the hot gas component. As shown in Fig. [5] 
the ejected gas mass exceeds the hot gas mass in low mass 
galaxies, and it is thus of crucial importance how its strip- 
ping is modelled. Here, we follow Font et al. (2008) and 
simply assume that in a given stripping event, the stripped 
fraction of ejected gas equals the stripped fraction of hot 
gas. If Rtrunc < /3Rvir, we assume that all the ejected gas is 
stripped. 

We now implement the above model into a SAM based 
on DLB07, which includes the first three modifications in 
section 5.2. No tidal stripping is included, and the AGN 
feedback efficiency is set to be the same as in DLB07 for 
easier comparison. We compare two different models: Rl) 
and R2). In model Rl), we use /3= 0.1, and in model R2), 
/3 =0.01. Results for the stellar mass bins log(M//i _1 M ) 
= 10-10.25 (top panel) and log(iW '/h' 1 M ) = 10.5-10.75 
(bottom panel) are shown as green dashed lines for model 
Rl) and as magenta dot-dashed lines for model R2) in Fig 
1141 For comparison, we also show results from DLB07 and 
our model 1), where the diffuse gas is stripped at the same 
rate as the dark matter. 

For low mass galaxies, we find that the ram-pressure is 
much more efficient at stripping out the gas than model 1), 
and both models Rl) and R2) produce passive fractions that 
are too high to match observations^ On the other hand, 

4 We have also tested a model with /9=0, which means that the 
central density of the hot gas is never affected by stripping. In this 
case, satellite galaxies have as high star formation rates as central 
galaxies. Interestingly, this indicates that it is not the reduced hot 
gas mass, but the reduced central density of the hot gas, which 
causes environmental effects in all our models. 
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Figure 14. Passive fraction as a function of cluster-centric radius 
at stellar masses of log(Af//i~ 1 MQ)=10 - 10.25 (top panel) and 
10.5 - 10.75 (bottom panel). Different linestyles and colours cor- 
respond to DLB07 (blue dotted lines), our model 1) (cyan solid 
lines) and the two models Rl) and R2) including ram-pressure 
stripping and no tidal stripping (green dashed lines and magenta 
dot-dashed lines) as indicated and described in the text. Crosses 
with errorbars show the observational results. Only clusters with 
masses above 10 14 /i -1 Af Q are used. 



in high mass galaxies, the ram-pressure effects are weaker 
than the tidal effects, and the ram-pressure models cannot 
produce a strong enough change in passive fraction as a 
function of clustercentric radius. In Fig. 1151 we show the 
distribution of SSFR for the stellar mass bin log(M/h~ 1 MQ) 
= 10-10.25. In magenta, we show results for the model R2). 
Clearly, there are too many galaxies with intermediate SSFR 
in this model. 

In summary, we find three problems with a model in 
which environmental effects arise solely due to ram-pressure 
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Figure 15. The distribution of satellite SSFR at stellar masses 
of log(M/h~ 1 Mq)=10 - 10.25. Different colours correspond to 
DLB07 (blue lines), our model 1) (cyan lines) and the ram- 
pressure stripping model R2) as indicated and described in the 
text. Black lines show the observational results. Only clusters with 
masses above 10 14 /i -1 M s are used. 



stripping. First, stripping effects in clusters for low mass 
galaxies are too strong to explain the observations. Second, 
stripping effects are not strong enough to reproduce the ob- 
servations for high mass galaxies. Third, there are too many 
low mass galaxies with intermediate SSFR. The first prob- 
lem could indicate that our understanding of the so-called 
"ejected phase" , which is important for low mass galaxies, is 
inaccurate. Although the ejected gas is currently unavailable 
for star formation, it maybe nevertheless strongly bound to 
the galaxy, and difficult to strip. Another potential reason 
for the problem is that SAMs like DLB07 are likely to overes- 
timate the hot gas mass in groups with masses between 10 12 
and lO 14 /i -1 M0, especially in the central, most dense re- 
gions of these systems (Bower et al. 2008) . This could lead to 
too strong pre-processing of cluster galaxies in groups. The 
second problem could indicate that tidal effects are in fact 
important for higher mass galaxies. By decreasing the dark 
matter content in subhaloes, tidal effects make high mass 
galaxies more vulnerable to ram-pressure stripping, but also 
to AGN feedback. The third problem indicates that ram- 
pressure stripping as implemented here does not produce 
a bimodal satellite galaxy population. A similar problem 
seems to be present in the ram-pressure stripping model of 
Font et al. (2008), as pointed out by Balogh et al. (2009). 
This could also hint at the importance of tidal effects, which 
we have shown to produce bimodal distribution in the SSFR 
of galaxies. 

We note that Font et al. (2008) have found that their 
ram-pressure model can reproduce the colours of satellite 
galaxies well. As they do not follow the orbits of individual 
satellites in detail, and use different assumptions for the hot 
gas profiles and the treatment of cooling in satellites, it is 
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hard to directly compare the two models and to pinpoint 
the reason for this discrepancy. 



8 SUMMARY 

We explore a set of modifications to the SAM of DLB07 with 
the goal to reproduce (i) the passive fractions as a function 
of cluster-centric radius out to the field, (ii) the detailed 
distribution of the SSFRs of satellite galaxies and (iii) the 
passive fractions of central galaxies as a function of stellar 
mass. We apply a realistic cluster finder to the SAM to allow 
a fair comparison with observational results from the cluster 
catalogue of vdL07. Our main results are the following. 

• A model in which the diffuse gas is stripped at the same 
rate that the dark matter subhalos lose mass due to tidal 
stripping can reproduce observations reasonably well. "Star- 
vation" is usually defined as the fast removal of the hot gas 
halo surrounding satellite galaxies (Larson, Tinsley & Cald- 
well 1980) due to both tidal and ram-pressure stripping. In 
this work, we find indications that this picture is in need 
of revision. First, stripping of the diffuse gas reservoir is 
slow. Second, it is crucial for the success of our model that 
even the gas ejected by supernovae is not stripped instanta- 
neously, but viewed as part of the overall gas reservoir. 

• After infall, the median subhalo identified at z=0 loses 
15-20% of its initial dark matter halo mass per Gyr, which 
implies that tidal stripping is typically only complete after 
5-7 Gyrs. Dark matter stripping seems to proceed linearly, 
and not exponentially. This is crucial for reproducing the 
observed bimodality in the specific star formation rates of 
satellites in our model. We also find that stripping of dark 
matter becomes more efficient at high redshift. 

• We present a simple method for implementing gas strip- 
ping in SAMs which do not follow the dark matter subhaloes 
of satellites. Our prescription is that 10-20% of the hot and 
ejected gas recorded at infall should be stripped per Gyr. 
This model gives good agreement with observations at z=0. 

• A model which only includes a physically motivated 
prescription for ram-pressure stripping of the hot gas in 
satellite galaxies seems to lead to an overproduction of low 
mass passive galaxies in clusters, but also too little environ- 
mental effects for high mass galaxies. 

• We have found that DLB07 do not reproduce the pas- 
sive fraction of central galaxies as a function of stellar mass 
in detail. In particular, DLB07 predict too few passive cen- 
tral galaxies with low masses, and too many with high 
masses. We find that we can improve agreement with ob- 
servations by reducing the AGN feedback efficiency, and by 
assuming that satellite galaxies in host haloes with masses 
below 1O 12 /i _1 M0 retain all their hot and ejected gas. How- 
ever, the reduction of AGN feedback leads to an overproduc- 
tion of the total number of massive and luminous galaxies. 
In order to simultaneously reproduce the low abundance, 
and the relatively high active fraction of massive galaxies, 
the form of the AGN feedback in DLB07 may need to be 
changed. One possible way to solve this problem might be to 
introduce some stochasticity into the shut-down of cooling at 
high halo masses, for example by taking into account the ef- 
fect of major mergers or morphological quenching (Hopkins 
et al. 2008; Martig et al. 2009). We also find that reproducing 
both the red luminosity function and the passive fraction of 



central and satellite galaxies is surprisingly difficult, indicat- 
ing either an inconsistency between different observations, 
or a potential problem with the dust treatment in the SAM. 
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